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Abstract:  The  possibility  for  controlling  both  the  probe-field  optical 

gain  and  absorption,  as  well  as  photon  conversion  by  a  surface-plasmon- 
polariton  near  field  is  explored  for  a  quantum  dot  located  above  a  metal 
surface.  In  contrast  to  the  linear  response  in  the  weak-coupling  regime, 
the  calculated  spectra  show  an  induced  optical  gain  and  a  triply-split 
spontaneous  emission  peak  resulting  from  the  interference  between  the 
surface-plasmon  field  and  the  probe  or  self-emitted  light  field  in  such  a 
strongly-coupled  nonlinear  system.  Our  result  on  the  control  of  the  mediated 
photon-photon  interaction,  very  similar  to  the  ‘gate’  control  in  an  optical 
transistor,  may  be  experimentally  observable  and  applied  to  ultra-fast 
intrachip/interchip  optical  interconnects,  improvement  in  the  performance 
of  fiber-optic  communication  networks,  and  developments  of  optical  digital 
computers  and  quantum  communications. 
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1.  Introduction 

It  is  well  known  that  photons  inherently  do  not  interact  with  each  other.  In  classical  electro¬ 
dynamics,  Maxwell  equations  are  linear  and  cannot  describe  any  photon-photon  interaction. 
However,  effective  photon-photon  coupling  could  exist  in  a  mediated  way,  e.g.,  through  their  in¬ 
teractions  with  matter.  Very  recently,  an  experiment  [1],  which  involves  firing  pairs  of  photons 
through  an  ultra-cold  atomic  gas,  provided  the  evidence  for  an  attractive  interaction  between 
photons  to  form  so-called  ‘molecules’  of  light.  In  general,  if  the  interaction  between  photons 
and  matter  is  strong,  the  optical  response  of  matter  will  become  nonlinear  and  the  resulting 
bandedge  optical  nonlinearities  [2]  will  enable  an  effective  photon-photon  interaction  [3].  An 
optical  transistor  [4]  could  be  based  on  this  idea,  where  ‘gate’  photons  control  the  intensity  of 
a  light  beam  ‘source’.  Optical  transistors  could  be  applied  to  speed  up  and  improve  the  perfor¬ 
mance  of  fiber-optic  communication  networks.  Here,  all-optical  digital  signal  processing  and 
routing  is  fulfilled  by  arranging  optical  transistors  in  photonic  integrated  circuits,  and  the  signal 
loss  during  the  propagation  could  be  compensated  by  inserting  new  types  of  optical  amplifiers. 
Moreover,  optical  transistors  are  expected  to  play  an  important  role  in  the  development  of  an 
optical  digital  computer  or  quantum-encrypted  communication. 

Most  previous  research  on  optical  properties  of  materials,  including  optical  absorption,  in¬ 
elastic  light  scattering,  and  spontaneous  emission,  used  a  weak  probe  field  as  a  perturbation 
to  the  studied  system  [5].  In  this  weak-coupling  regime,  the  optical  response  of  electrons  de¬ 
pends  only  on  the  material  characteristics  [6],  and,  therefore,  no  photon-photon  interaction  is 
expected.  However,  the  strong-coupling  regime  could  be  reached  with  help  from  microcavi¬ 
ties,  and  the  experimental  effort  on  searching  for  polariton  condensation  (resulting  from  strong 
light-electron  interaction)  in  semiconductors  continues  to  produce  results  [7-9].  A  general  re¬ 
view  of  exciton-polariton  condensation  may  be  found  in  [10],  Successful  demonstration  of 
room-temperature  polariton  lasing  without  population  inversion  in  semiconductor  microcavi¬ 
ties  using  both  optical  pumping  [1 1, 12]  and  electrical  injection  [13, 14]  have  made  it  possible 
to  achieve  ultra-low  lasing  thresholds  and  very  small  emitter  sizes  comparable  to  the  emitted 
wavelength.  Semiconductor  exciton-polariton  nanolasers  could  advance  intrachip  and  interchip 
optical  interconnects  by  integrating  them  into  semiconductor-based  photonic  chips.  They  might 
also  have  applications  in  medical  devices  and  treatments,  such  as  spatially  selective  illumina¬ 
tion  of  individual  neuron  cells  to  locally  control  neuron  firing  activities  in  optogenetics  and 
neuroscience,  and  near-field  high-resolution  imaging  beyond  the  optical  diffraction  limit  as 
well. 

Theoretically,  a  big  hurdle  also  exists  for  studying  photon-photon  interactions  in  the  strong¬ 
coupling  regime  mainly  due  to  intractable  numerical  calculations  for  systems  with  very  strong 
nonlinearity.  The  existence  of  strong  nonlinearity  in  such  a  system  means  that  any  perturbative 
theories,  e.g.  using  bare  electron  states  or  linear  response  theory  [5],  become  inadequate  for 
describing  both  field  and  electron  dynamics  in  this  system.  The  presence  of  an  induced  polar¬ 
ization,  regarded  as  a  source  term  in  Maxwell  equations  [15, 16],  from  photo-excited  electrons 
makes  it  impossible  for  us  to  solve  the  field  equations  by  simply  using  finite-element  anal¬ 
ysis  [17]  or  finite -difference-time-domain  methods  [18].  Although  the  semiconductor-Bloch 
equations  [19]  and  density-matrix  equations  [5,20],  derived  from  many-body  theory,  are  able 
to  accurately  capture  the  nonlinear  optical  response  of  electrons,  the  inclusion  of  pair  scatte- 
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ring  effects  on  both  energy  relaxation  and  optical  dephasing  precludes  an  analytic  approach  for 
seeking  a  solution  of  these  equations.  Consequently,  there  exist  only  very  few  theoretical  studies 
on  resonantly-coupled  excitons  to  microcavity  field  [21],  which  depend  heavily  on  computer 
simulations,  and  focus  on  simplified  one-dimensional  strongly-coupled  microcavity  systems. 

Physically,  not  only  high-quality  microcavities  [22]  but  also  intense  surface-plasmon  near 
fields  [23, 24]  should  be  employed  for  reaching  the  strong-coupling  goal  in  semiconductors. 
In  this  paper,  we  solve  the  self-consistent  equations  for  strongly-coupled  electromagnetic -field 
dynamics  and  electron  quantum  kinetics  in  a  quantum  dot  above  the  surface  of  a  thick  metal¬ 
lic  film.  This  is  a  situation  that  has  not  been  fully  explored  so  far  from  either  a  theoretical  or 
experimental  point  of  view.  This  is  done  based  on  finding  an  analytical  solution  to  Green’s  func¬ 
tion  [25,26]  for  a  quantum  dot  coupled  to  a  semi-infinite  metallic  material  system,  which  makes 
it  easy  to  calculate  the  effect  of  the  induced  polarization  field  as  a  source  term  in  Maxwell  equa¬ 
tions.  In  our  formalism,  the  strong  light-electron  interaction  is  reflected  in  the  photon-dressed 
electronic  states  with  a  Rabi  gap  and  in  the  feedback  from  the  induced  optical  polarization  of 
dressed  electrons  to  the  incident  light.  The  formalism  derived  in  this  paper  goes  beyond  the 
weak-coupling  limit  and  deals  with  a  much  more  realistic  structure  in  the  strong-coupling  limit 
for  the  development  of  a  surface-plasmon  polariton  laser  with  a  very  low  threshold  pumping. 
Our  results  clearly  demonstrate  the  ability  to  control  probe-field  optical  gain  and  absorption, 
and  photon  conversion  by  a  surface-plasmon  field  with  temperature -driven  frequency  detuning 
in  such  a  nonlinear  system  led  by  dressed  electron  states,  very  similar  to  the  ‘gate’  control  in  an 
optical  transistor.  These  conclusions  should  be  experimentally  observable  [27,28].  On  the  other 
hand,  our  numerical  results  also  provide  an  example  that  demonstrates  the  so-called  quantum 
plasmonics  [29],  where  the  properties  of  surface-plasmon  polaritons  and  quantum-confined 
electrons  are  hybridized  through  near-field  coupling. 

In  Sec.  2,  we  introduce  our  physics  model  and  derive  self-consistent  equations  for  determin¬ 
ing  the  coupled  scattering  dynamics  of  a  surface-plasmon  field  and  the  quantum  kinetics  of 
electrons  in  quantum  dots.  Section  3  is  devoted  to  a  full  discussion  of  our  numerical  results, 
including  scattering  and  optical  absorption  of  the  surface-plasmon-polariton  field  by  quantum 
dots,  spontaneous  emission,  and  the  nonlinear  optical  response  of  the  dressed  electron  states. 
Some  concluding  remarks  are  given  in  Sec.  4. 

2.  Model  and  theory 

Our  model  system,  as  shown  in  Fig.  1,  consists  of  a  semi-infinite  metallic  material  with  a  semi¬ 
conductor  quantum  dot  above  its  surface.  A  surface-plasmon-polariton  (SPP)  field  is  locally 
excited  through  a  surface  grating  by  normally-incident  light.  This  propagating  SPP  field  fur¬ 
ther  excites  an  interband  electron-hole  (e-h)  plasma  in  the  quantum  dot.  The  induced  optical- 
polarization  field  of  the  photo-excited  e-h  plasma  is  strongly  coupled  to  the  SPP  field  to  produce 
split  degenerate  e-h  plasma  and  SPP  modes  with  an  anticrossing  gap.  A  brief  description  of  our 
self-consistent  formalism  was  reported  earlier  [16].  In  order  to  let  readers  follow  easily  the 
details  of  our  model  and  formalism,  we  present  here  the  full  derivation  of  the  Maxwell-Bloch 
numerical  approach  for  an  SPP  field  coupled  to  a  photo-excited  e-h  plasma  in  the  quantum  dot. 

2.1.  General  formalism 

The  Maxwell  equation  for  a  semi-infinite  non-magnetic  medium  in  position-frequency  space 
can  be  written  as  [25] 


e0c2 


^loc(r;  to)  , 


(1) 
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Fig.  1.  Schematic  illustration  for  a  semi-infinite  metal  and  a  quantum  dot  above  its  surface 
atx3  =  0.  Flere,  the  surface-plasmon  polariton  (SPP)  is  locally  excited  by  incident  light  with 
the  help  of  a  surface  grating.  The  propagating  SPP  field  further  excites  e-h  pairs  (plasmas) 
in  the  adjacent  quantum  dot.  As  a  result,  the  optical-polarization  field  of  the  photo-excited 
e-h  plasma  is  strongly  coupled  to  the  propagating  SPP  field  to  form  split  plasma-SPP  modes 
with  an  anticrossing  gap.  Also,  a  probe-field  is  used  for  studying  the  photon  dressing  effect. 


where  E(r;co)  is  the  electric  component  of  an  electromagnetic  field, 
H{ r;  co)  =  —  ( -  ]  Vx£(r;  co)  is  the  magnetic  component  of  the  electromagnetic  field, 

V  &>At()  / 

r  =  (xi ,  X2  .  x3 )  is  a  three-dimensional  position  vector,  co  is  the  angular  frequency  of  the  incident 
light,  £o,  fio  and  c  are  the  permittivity,  permeability  and  speed  of  light  in  vacuum,  ,S^loc('r;  co)  is 
an  off-surface  local  polarization  field  generated  by  optical  transitions  of  electrons  in  a  quantum 
dot,  and  the  position-dependent  dielectric  function  is 


£b(x3;  co)  = 


£d  ,  for  x'3  >  0 
£m(co)  ,  for  x'3  <0 


(2) 


Here,  £<i  characterizes  for  the  semi-infinite  dielectric  material  in  the  region  x3  >  0,  while  £m(&>) 
is  the  dielectric  function  of  the  semi-infinite  metallic  material  in  the  region  X3  <  0.  For  the 
Maxwell  equation  in  Eq.  (1),  we  introduce  the  Green’s  function  Sffiv(r,r';  co)  that  satisfies  the 
following  equation 


£b(x3;  co 


CO- 


- 


dx^dxn 


(r,r';  co)  =  8XvS(r-r') 


(3) 


where  ^  0  is  the  Laplace  operator,  <5;  u  represents  the  Kronecker  delta,  and  the  indices 

A,  jtt  =  1, 2,  3  indicate  three  spatial  directions.  Using  the  Green’s  function  defined  in  Eq.  (3), 
we  can  convert  the  Maxwell  equation  in  Eq.  (1)  into  a  three-dimensional  integral  equation 


2 

Eu(r;  co)  =E^u)(r;  co )--^£  / d3r'%v(r,r'-,  co)  ^“(r';  co)  , 
£0C  y  J 


where  E^>:  (r;  co)  is  a  solution  of  the  corresponding  homogeneous  equation 


eb(x3,  co)  8^v  ^  +5^vVr 

/i  V 


40)(r;  <o)  =  0, 


(4) 


(5) 


and  the  source  term  ,'3'a|/0C(V;  co)  generally  depends  on  the  electric  field  in  a  nonlinear  way  and 
can  be  determined  by  the  Bloch  equation  [2, 15]. 
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2.2.  Solving  for  Green ’s  function 

For  a  semi-infinite  medium,  the  Green’s  function  can  be  formally  expressed  by  its  Fourier 
transform 


% 


gv 


(r,r';  to)  =  J 


d2k  || 

(2k)2 


'kn-Crirrii) 


^v(k||,0)|x3,X3)  , 


(6) 


where  we  have  introduced  the  two-dimensional  vectors  rii  =  (xi,X2)  and  k||  =  Substi¬ 

tuting  Eq.  (6)  into  Eq.  (3),  we  obtain 


or  ,  , 

£b^~kl  +  dx2 


k\ki 


— ik\ 


dx  3 


k\k2  -iki~— 

dx  3 

(O2  t  d2  d 

£b  — —  k\  +  ~z~2  —ikl  -7 — 

c-  dx 3  dx  3 

d  CO2  2 

lk,2  t  c 2 


'  dx  3 
=  5(x3  -x'3) 


gll  g  12  gl3 

821  g22  823 

8 31  8 32  t?33 


1  0  0 
0  1  0 
0  0  1 


After  a  rotational  transformation  [25]  is  performed  in  k||  -space,  i.e., 

//iv(*||.®  1x3,4)  =  S  ^'(k!|)^vv'(k|j)^v(k||,co|x3,X3)  , 


n  y 

where  the  rotational  matrix  «^(kn)  is 

^(kn)  =  4 

we  acquire  an  equivalent  version  of  Eq.  (7) 

co2 


k  i  &2  0 

— k2  k\  0 

0  0  /til 


£b 


•'2  t/2 


-  + 


0 


dx2 


0 


—iku 


— iku 


dxi, 


of  ,2  d2 

£b^~kU  + 

0 


=  5(x3  -x'3) 


<iX3 


£0  ,2 


1  0  0 
0  1  0 
0  0  1 


fn 

/12 

/l3 

/21 

fl2 

fl3 

/31 

hi 

/33 

(7) 


(8) 


(9) 


(10) 


To  obtain  the  solution  of  Eq.  (10),  we  need  to  employ  both  the  finite- value  boundary  condition 
at  X3  =  ±00  and  the  continuity  boundary  condition  at  the  X3  =  0  interface.  This  leads  to  the 
following  five  non-zero  /juV(A||,fi>|x3,X3)  elements  [25,26]: 


/22(*||>®|x3,x'3) 
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i  \  2  p 


2pJ  Pd+P 


eiPA*3-ipx3 


1  '  M*3-A I  _  .-<>te+4) 


= 


2Pj  L  Pd  +  P 

,i>dl-«3-4l  -|_  l?d~  P  e'7>d(*3+-4) 

2pdy  r  Pd+p 


_  ( J_\  2Pd  e-ipte-^) 
\2pd  /  Pd  +  P 


X3  >  0  ,  X3  <  0 
,  X3  <  0  ,  X3  <  0 
,  X3  >  0  ,  X3  >  0 
X3  <  0  ,  X3  >  0 


(11) 


/l3(fc||,  CO |x3,X3) 

ik\\c2  r  2  eM(®)pd 
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2eM(ft))ffl2  LeM(®)pd  +  £dP. 

ik\,c2 


ipix-i-ipx^ 


2  eM(ffl)ffl2 

ifeii  c2 


‘>1-3-41  sgn(x3  —  X3)  +  £^m\PA~edP  e-ip{x 3+4)' 

eM(co)/3d  +  edP 


2edco2 

//hie2 


Odl-3-41  sgn(x3  -y,)  +  e‘>d(-3+4)' 

£M(®)pd  +  £dP 


2edp 


2ed ft)2  L£M(a>)pd  +  £dP_ 


-i/M3+ipd4 
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(12) 
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f33{k\\,CO\x3,x'3) 
ik},c2  r 
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JPd^-ipA 
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^  (x3  X3  ) 
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,'>1-3-41  _  £M(ft>)Pd~£dP  e-ip(x 3+y3) 
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2^2 


£d  ft) 
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2pd£d  co2 


,'>dl— 3— 41  -4-  £MMPd~£dP  jPd(^+4) 


£M(®)/7d  +  £dP 


CO2  L£M(®)Pd  +  £dP. 


-/pjr3+ipd4 
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X3  <  0  ,  X3  <  0 

(13) 

X3  >  0  ,  X3  >  0 

X3  <  0  ,  x'3  >  0 


/n(^l|,  CO  |x3,x'3) 
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X3  >  0  ,  X3  <  0 


ipdpc  [ _ 1 

ft)2  [eM(ft))f>d  +  ed^. 


ipiX^-ip^ 


ipc~ 


2  £M{(0)(D2 


,ip\x3-x!3\  £M(<0)Pd-£dP  e-ip(x 3+J/3) 

£M(®)/?d  +  £dP 
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ipdc- 
2  £dft>2 

ip<jC2 

2  edco2 


Bfate-4I  _  £M(^)Pd-£dP 


2£dp 


£M(ft>)/?d  +  £dP_ 


£M(ft))Pd  +  £dP 


-ipxi+ipA^ 


X3  <  0  ,  X3  <  0 
X3  >  0  ,  X3  >  0 
X3  <  0  ,  X3  >  0 


(14) 


/3l(£j|,ft)|x3,X3) 

/Arn  c2 


ft)* 


£M(®)pd  +  £dP. 


ipix3-ipxl3 


ihC  f  >‘>1*3-41  sgn(x3  -  x' )  -  e-‘>(*3+4)' 

£]Vl(ft))Pd  +  £dP 
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2eM(ffl)®2 

^  r  fate-41  /  _ y  )  _  £M(to)Pd-£dP  efate+V3) 

eM(®)pd  +  £dP 


2edco2 

;liic2  r 


ft)* 


Pd 


£M(®)Pd  +  £dP. 


-ipx3+ipix'3 


where  sgn(x)  is  the  sign  function. 


Pd(k\\,(o)  =  \/ £d^  -kj 


X3  >  0  ,  X3  <  0 
X3  <  0  ,  X3  <  0 

(15) 

X3  >  0  ,  X3  >  0 
X3  <  0  ,  X3  >  0 


(16) 


P(k\U°>)  =  y  £m{(0)  ^2  -  kj  ,  (17) 

Im[pd(^|| ,  ft))]  >  0  and  lm\p{k\\ ,  ft))]  >  0.  In  addition,  from  these  non-zero  fpV(k\\ ,  (O^Xj)  func¬ 
tions,  we  obtain 


^v(k||,ft)|x3,x,3)  =  £  /Mv(^||,ft)|x3,X3)^/M(k||)^v/v(k||)  ,  (18) 

n'y 

which  can  be  substituted  into  Eq.  (6)  to  calculate  the  Green’s  function  (fuv(r.  r‘;  co)  in  position 
space. 

2.3.  Local  polarization  field 

In  order  to  find  the  explicit  field  dependence  in  £Ploc(r;  ft)),  we  now  turn  to  the  study  of 
electron  dynamics  in  a  quantum  dot.  Here,  the  optical-polarization  field  .^loc(r;  co)  plays 
a  unique  role  on  bridging  the  classical  Maxwell  equations  for  electromagnetic  fields  to  the 
quantum-mechanical  Schrodinger  equation  for  electrons.  The  electron  dynamics  in  photo- 
excited  quantum  dots  can  be  described  quantitatively  by  the  so-called  semiconductor  Bloch 
equations  [30-32].  These  generalize  the  well-known  optical  Bloch  equations  in  two  aspects 
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including  the  incorporation  of  electron  scattering  with  impurities,  phonons  and  other  electrons 
as  well  as  many-body  effects  on  dephasing  in  the  photo-induced  optical  coherence. 

The  physical  system  considered  in  this  paper  is  illustrated  in  Fig.  2,  where  we  assume  two 
levels  for  electrons  and  holes,  respectively,  in  a  quantum  dot.  These  two  energy  levels  of  both 
electrons  and  holes  are  efficiently  coupled  by  phonon  scattering  at  high  temperatures.  Addi¬ 
tionally,  the  lowest  electron  and  hole  energy  levels  are  optically  coupled  to  each  other  by  an 
incident  SPP  field  to  form  the  dressed  states  of  excitons.  The  SPP-controlled  optical  properties 
of  quantum-dot  excitons  can  either  probed  by  a  plane-wave  field  or  seen  from  the  spontaneous 
emission  of  excitons. 


"2 

»e 
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...■■■"'phonon  coupling 


G 

EP 
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...°L 

spontaneous 
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®em 

l  zz 
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evanescent 
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Fig.  2.  Schematic  illustration  for  a  system  incorporating  the  generation  of  quantum-dot 
excitons  by  a  SPP  field  with  frequency  msp  and  probed  by  a  plane-wave  field  Ep  with 
frequency  cop.  Here,  Sq  is  the  energy  bandgap  of  the  host  semiconductor,  ef  and  e*l  stand 
for  the  energy  levels  of  electrons  and  holes,  respectively,  with  i ,  j  =  1, 2,  •  ■  •.  In  addition, 
(Oem  represents  the  frequency  of  spontaneous  emission  from  photo-excited  excitons,  and 
the  ground  states  of  electrons  and  holes  are  coupled  to  their  first  excited  states  by  lattice 
phonons  at  finite  temperatures. 


For  photo-excited  spin-degenerated  electrons  in  the  conduction  band,  the  semiconductor 
Bloch  equations  with  t  =  1 , 2,  •  •  •  are  given  by 


dn  e  2 

— t  =  -  >  Im 
dt  h  ■ 


dne 

dt 


A: .  l  -^sp i  n  i  ■ 


(19) 


rel 


where  t%sv  is  the  spontaneous  emission  rate  and  if  represents  the  electron  level  population.  In 


Eq.  (19),  the  term  marked  ‘rel’  is  the  non-radiative  energy  relaxation  for  rf(,  and  the  Yjl ,  •dtf'j, 
and  Vfj.j  f  terms  are  given  later  in  the  text. 

Similarly,  for  spin-degenerate  holes  in  the  valence  band,  the  semiconductor  Bloch  equations 
with  j  =  1 , 2,  •  •  •  are  found  to  be 


A 

dt 


XIm[(T/)*(,< 


A 

dt 


—  <5/.1  &sp n \n\  i 

rel 


(20) 


where  nh-  stands  for  the  hole  energy  level  population.  Again,  the  non-radiative  energy  relaxation 
for  ifj  is  incorporated  in  Eq.  (20).  Moreover,  we  know  from  Eqs.  (19)  and  (20)  that 
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(21) 


Ne(t)=2^nee(t)=2^n)(t)=Nh(t), 

t  j 

where  Ne(t)  and  N^(t)  are  the  total  number  of  photo-excited  electrons  and  holes,  respectively, 
in  the  quantum  dot  at  time  t. 

Finally,  for  spin-averaged  e-h  plasmas,  the  induced  interband  optical  coherence,  which  is 
introduced  in  Eqs.  (19)  and  (20),  with  j  =  1,  2,  ■  •  •  and  /  =  1 . 2 ,  ■  ■  ■  satisfies  the  following  equa¬ 
tions. 


ih  —Yj  = 
dt  1 


£^(co)  +  £j(co)  —  h((o  +  iyo)  Y^+(l—  nf  — 

V  „h  ( T/hh  T/hh  \  V  T/eh 

yjJi'JiJ  ^jnh^hJJA 


+Y 


+Y 


X  nh  (Vhh;h/  -  vhhAh )  "  X  n)i  vth\h / 
M  h 


(22) 


where  tijo  =  hy^  -\-  tiyext  is  the  total  energy-level  broadening  due  to  both  the  finite  carrier  life¬ 
time  and  the  loss  of  an  external  evanescent  field,  CO  is  the  frequency  of  the  external  field,  and 
e^(co)  and  £j(co)  are  the  kinetic  energies  of  dressed  single  electrons  and  holes,  respectively  (see 
Appendix  A  with  a  =  1).  In  Eq.  (22),  the  diagonal  dephasing  (yo)  of  Yl,  the  renormalization 
of  interband  Rabi  coupling  (K/  vfj.  ■  e),  the  renormalization  of  electron  and  hole  energies  (third 
and  fourth  terms  on  the  right-hand  side),  as  well  as  the  exciton  binding  energy,  are  all  taken 
into  consideration.  Since  the  e-h  plasmas  are  independent  of  spin  index  in  this  case,  they  can  be 
excited  by  both  left-circularly  and  right-circularly  polarized  light.  The  off-diagonal  dephasing 
of  y/  has  been  neglected  due  to  low  carrier  density. 

The  steady-state  solution  to  Eq.  (22),  i.e.  under  the  condition  of  dY](  /dt  =  0,  is  found  to  be 


Yl(t\m) 


1  -n\{t)-n)(t) 

— eh 

h(co  +  iyo)  —  TiQ.t  ;-(co|f) 


(23) 


where  the  photon  and  Coulomb  renormalized  interband  energy-level  separation  hSl^j(co\t)  is 
given  by 


hUr;j(co\t)  =  £e({co\t)+£){(o\t)-Veenj.jti 


reh 


h 


h  hi*  hii 

The  steady-state  solution  in  Eq.  (23)  can  be  substituted  into  Eqs.  (19)  and  (20)  above. 

The  Coulomb  interaction  matrix  elements  introduced  in  Eqs. (19),  (20)  and  (22)  are  defined 
as 

VhAAA=J  d\jd\'  [v£(r)]*  [v^OO]*  4^|r_r/|  Vt3  OO  (r) 


87T2£o£i 
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Vhj2j3,j4=  j  ^  j  ^ 


Sn2£o£b 


Id2q> 


’*<■'>]  SSSF^TitOTiM 
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\  JUJ2'J3,J4  J 


(26) 


vtrJ'/'  =  jd^rj dW'  [^(r)]' 


8n2£o£b 


1 


q+q.s 


(27) 


where  the  static  screening  length  1  /qs  at  temperatures  (k/jT  2>  £»  is  determined  from 


(28) 

Here,  5?  is  the  cross-sectional  area  of  a  quantum  dot,  T  is  the  lattice  temperature,  i//f(r)  and 
t j/-  (r)  are  the  envelope  wave-functions  of  electrons  and  holes  in  a  quantum  dot  (see  Appendix 
A),  and  £/,  is  the  average  dielectric  constant  of  the  host  semiconductor.  The  two  dimension¬ 
less  form  factors  (see  Appendix  A)  introduced  in  Eqs.  (25)-(27)  for  electrons  and  holes  due  to 
quantum  confinement  by  a  quantum  dot  are  defined  by 


^>(q||)  =e-qW^j d\  [vff(r||)]*eiqll’rH  ^(r„) 


(29) 


(30) 


where  Jzfo  is  the  thickness  of  a  disk-like  quantum  dot.  In  addition,  the  matrix  elements  em¬ 
ployed  in  Eqs.  (19),  (20)  and  (22)  for  the  Rabi  coupling  between  photo-excited  carriers  and  an 
evanescent  external  field  E(r;t)  =  0(t)E (r;  ft))  e~,mt  are  given  by 


-8(a  Sj.  i  G(t) 


Efj(co)-  dCjV 


(31) 


where  0  (x)  is  a  unit  step  function,  the  static  interband  dipole  moment  dc,v  (see  Appendix  A)  is 


dc,v  =  J  d3 r  [uc  (r)]  *  r  «v  (r)  =  d*  v  ,  (32) 

mc( r)  and  uv(r)  are  the  Bloch  functions  associated  with  conduction  and  valence  bands  at  the 
T-point  in  the  first  Brillouin  zone  of  the  host  semiconductor,  and  the  effective  electric  field 
coupled  to  the  quantum  dot  is 


=  J d3r  [i/r|(r)]*E(r;  c)  [$(r) 


(33) 


The  Boltzmann-type  scattering  term  [33]  for  non-radiative  electron  energy  relaxation  in 
Eq.  (19)  is 


dnt 

dt 


rel 


^(in)(l_ne)_^(out)ne! 


where  the  microscopic  scattering-in  and  scattering -out  rates  are  calculated  as 


(34) 
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(36) 


Here,the  primed  summations  in  Eqs.  (35)  and  (36)  exclude  the  terms  satisfying  either  j  =  f 
or  l'  =  l,  tylph(Qo)  =  [exp (ftQo/^B r)  —  l]~x  is  the  Bose  function  for  the  thermal-equilibrium 
phonons,  and  fio  and  Tph  are  the  frequency  and  lifetime  of  longitudinal-optical  phonons  in 
the  host  semiconductor.  Similarly,  the  Boltzmann-type  scattering  term  for  hole  non-radiative 
energy  relaxation  in  Eq.  (20)  is 


dn^ 

~dt 


J  =Wf\\-n))-W^]n) 

rel 


where  the  scattering-in  and  scattering-out  rates  are 
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and  again  the  primed  summations  in  Eqs.  (38)  and  (39)  exclude  the  terms  satisfying  either 
/  =  j  or  £  =  (! .  The  coupling  between  the  longitudinal-optical  phonons  and  electrons  or  holes 
in  Eqs.  (35),  (36),  (38)  and  (39)  are 


veP 

vt,l< 


hp 


11' 


e2h£i.Q 

87T2£o 

e2hQ.Q 

8k2£q 


ry 


+qs 

l 

q\\  +q. s 


(40) 


(41) 


where  E,x  and  es  are  the  high-frequency  and  static  dielectric  constants  of  the  host  polar  semi¬ 
conductor. 

By  generalizing  the  Kubo-Martin-Schwinger  relation  [20],  the  time-dependent  spontaneous 
emission  rate,  ^sp(f),  introduced  in  Eqs.  (19)  and  (20),  can  be  expressed  as 
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(42) 


(43) 


£q(T)  =  <%q(0)  —  5.41  x  1(T4  T2/(T  +  204)  (in  units  of  eV)  is  the  energy  bandgap  of  the  host 
semiconductor,  po(co)  =  or /c3n2h  is  the  density-of-states  of  spontaneously-emitted  photons 
in  vacuum,  mo  is  the  free  electron  mass,  m*  is  the  effective  mass  of  electrons,  and  the  Coulomb 
renormalization  of  the  energy  bandgap  Sc(t)  is  found  to  be 

«(')  =  I »!,(')  M«„> ~vm)  +  £»;,(<)  (<+„,  -Vf5, ;1Jl) 


h 


iJiOid 


1  -n\{t)-n\{t) 


i/eh 

1:1,1 


(44) 


In  Eq.  (44),  the  first  two  terms  are  associated  with  the  Hartree-Fock  energies  of  electrons  and 
holes,  while  the  remaining  terms  are  related  to  the  exciton  binding  energy. 

Finally,  the  photo-induced  interband  optical  polarization  ,'+loc(r;  to),  which  is  related  to  the 
induced  interband  optical  coherence,  by  dressed  electrons  in  the  quantum  dot  is  given  by  [2] 


^loc(r;  co)  =  2  ||(r)|' 


dcv{/ 


d3  r>1e(r')V'i(r') 


x  -  lim 

n  t->oo 


i-»i(Q-»i(0 


(45) 


where  dcv  =  r/c,v e,i  represents  the  interband  dipole  moment  [see  Eq.  (32)],  ec|  is  the  unit  vector 
of  the  dipole  moment,  and  | (r)  |2  comes  from  the  confinement  of  a  quantum  dot. 
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2.4.  Self-consistent  field  equation 

Since  the  wavelength  of  the  incident  light  is  much  larger  than  the  size  of  a  quantum  dot,  we 
can  treat  the  quantum  dot,  which  is  excited  resonantly  by  the  incident  light,  as  a  point  dipole  at 
r  =  ro  =  (0,0, zo),  he.  we  can  assume  ^loc(r';  to)  =  ^loc(to)  <5(r'  —  ro)  in  Eq.  (4)  to  neglect 
its  geometry  effect.  This  greatly  simplifies  the  calculation  and  gives  rise  to 

2 

£ju( r;  (o)=E^\ r;  to)  -  £  %v(r,rQ-  to)^°c(to)  ,  (46) 

where 
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to  +  iyo-Qf^tolf) 


•<1(0 


-0(f)  [£(r0;  to)  -d, 


-]{/ 


^3ri^(r)vf(r) 


(47) 

(48) 


Substituting  Eqs.  (47)  and  (48)  into  Eq.  (46),  we  get  the  following  nonlinear  equations  for  the 
electromagnetic  field 


Eu{ r;  to)  =£,(I0)(r;  to) 


2  to2 
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h  t->°° 


l-»i(Q-»i(0 
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(49) 


where  the  quantum-dot  level  populations  n|(f)  and  njl(f)  depend  nonlinearly  on  E(ro;  to)  in  the 
strong-coupling  regime. 

If  the  electromagnetic  field  is  not  very  strong,  we  can  neglect  the  pumping  effect.  In  this 
linear-response  regime,  we  can  write  down  the  electron  and  hole  populations  in  a  thermal- 
equilibrium  state  [without  solving  Eqs.  (19)  and  (20)] 
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(50) 


(51) 


where  fo(x)  is  the  Fermi  function,  and  and  fi h  are  the  chemical  potentials  of  electrons  and 
holes,  respectively,  determined  by  Eq.  (21).  As  a  result  of  Eqs.  (50)  and  (51),  we  get  from 
Eq.  (49)  the  linearized  self-consistent  field  equation  at  r  =  ro 
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(53) 
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«dX^vi  (ro,ro;£o)e^ 

Vi 


where,  according  to  Eq.  (6),  we  have 


r  c/2kti 

^v(ro,ro;®)  =  y  snv (k|| i ffl|z0>zo)  • 


(54) 


The  solution  £(ro;  ft) )  of  the  linear-matrix  equation  in  Eq.  (52)  can  be  substituted  into  Eq.  (49) 
to  yield  the  spatial  distribution  of  the  electromagnetic  field  E( r;  to)  at  all  positions  other  than 
r  =  r0,  i.e., 


X  ®)4'}(ro;  ®) 
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r')V'f(r') 


l-/o(£l)-/o(^) 
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Vl 


In  order  to  find  the  coupled  e-h  plasma  and  plasmon  dispersion  relation  ft)  =  Qex-pi(k|| ),  we 
Fourier  transform  both  £(r;  ft))  and  E,0':i  (r;  ft))  in  Eq.  (46)  with  respect  to  rn .  This  leads  to 


Efi  (k|| ,  fi)|x3)  =  E^u> (kj| ,  ft)|x3)  -  ~~^2  S  ^v(k|| ,  o|*3,zo)  3^c((i>)  ■ 
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(56) 


After  setting  x3  =  zq  in  Eq.  (56),  we  get 
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(57) 


Here,  the  vanishing  of  the  determinant  of  the  coefficient  matrix  in  Eq.  (57)  determines  the  cou¬ 
pled  e-h  plasma  and  plasmon  dispersion  relation  ft)  =  £2ex-pi(k||).  We  emphasize  that  the  as¬ 
sumption  of  thermal-equilibrium  states  for  electrons  and  holes  is  just  for  obtaining  analytical 
expressions.  Therefore,  some  qualitative  conclusions  can  be  drawn  for  guidance  from  these  an¬ 
alytical  solutions.  Our  numerical  results,  however,  are  based  on  the  non-thermal-equilibrium 
steady  states  calculated  after  solving  self-consistently  the  coupled  Maxwell-Bloch  equations. 

By  assuming  an  incident  SPP  field  within  the  xixi-plane,  we  can  write 


£(0)(r;ft)sp)=£spe'k«(%P>D»^  [/koj33(Eo,ft)sp)-x3Eo(ft)sp)]  e^ko(^r)*ll  e-ft(*o,<%p)*3  ,  (58) 

where  xm  =  {xi ,  X2},  ko  and  x3  are  the  unit  vectors  in  the  ko  =  Eo(ft)sp){cos  Oq ,  sin 0o}  and  x3 
directions,  £sP  is  the  field  amplitude,  ft)sp  is  the  field  frequency,  Oq  is  the  angle  of  the  incident 
SPP  field  with  respect  to  the  xi  direction.  Do  =  {— xg,  —yg}  is  the  position  vector  of  the  surface 
grating,  and  the  two  wave  numbers  are 
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Pi  (*0,  ft^sp) 


(60) 


with  Re[^0(©sp)]  >  0  and  Rc\Ji-i,(k().  0 Js p) ]  >  0.  Here,  the  in-plane  wave  number  kb  is  produced 
by  the  surface-grating  diffraction  of  the  /^-polarized  normally-incident  light,  which  in  turn  de¬ 
termines  the  resonant  frequency  ft)  of  the  SPP  mode.  Equation  (59)  stands  for  the  full  dispersion 
relation  of  the  SPP  field,  including  both  radiative  and  non-radiative  parts.  From  Eq.  (58),  it  is 
easy  to  find  its  Fourier  transformed  expression 
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2.5.  Quantum-dot  absorption 

On  the  basis  of  the  electromagnetic  field  £(ro;  ft))  at  the  quantum  dot,  we  are  able  to  compute 
the  time-resolved  nonlinear  interband  absorption  coefficient  of  dressed  electrons  in  a  quantum 
dot  for  the  SPP  field  [34].  In  this  case,  we  find 
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where  aspp(ft);  t)  is  the  complex  Lorentz  function  given  by 
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and  the  scaled  refractive  index  function  nspp(ft);  t)  can  be  calculated  by 


(64) 


ftspp (CO,  t ) 
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In  Eqs.  (63)  and  (64),  the  dressed-state  effects  on  both  the  level  population  and  dipole  moment 
have  been  included.  In  addition,  we  have  introduced  the  following  notations  in  Eqs.  (63)  and 
(64) 
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2.6.  Probing  quantum-dot  dressed  states 

We  are  also  able  to  compute  the  time-resolved  linear  interband  absorption  coefficient  of  elec¬ 
trons,  dressed  by  the  SPP  field,  for  a  weak  probe  field  (not  the  strong  SPP  field)  on  the  basis 
of  the  electromagnetic  field  //(r(l;  co)  at  the  quantum  dot  calculated  above  [34].  Assuming  a 
spatially-uniform  probe  field  Ep{t)  =  0(t  —  z)Epe~,COpt  where  T  is  the  delay  time,  the  probe- 
field  absorption  coefficient  j3abs(&>p;  t)  of  the  lowest  dressed  state  is  given  by  Eq.  (62)  with  the 
replacements  of  ft),  )ispp,  and  aspp  by  cop,  npf,  and  Opf,  respectively,  where 
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Here,  using  Eq.  (24)  we  have 
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and 


Ti(»±(t)  =  /jft)±A(co;  0  •  (71) 

Moreover,  the  time-resolved  photoluminescence  spectrum  /7,em(ftj/;  t)  of  dressed  electrons 


in  the  quantum  dot  is  proportional  to 
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3.  Numerical  results  and  discussions 

3.7.  Results  for  the  dynamics  of  an  SPP  field 

In  the  first  part  of  our  numerical  calculations,  we  have  taken:  ,i/()  =  100  A,  L„  =  100  A,  ml  = 
0.067  mo,  mj(  =  0.62  mo,  Qq  =  45°,  =  yg  =  610  A,  =  12,  es  =  11,  £>,  =  13,  ft£2o  =  36meV 

and  7zyeh  =  =  ftyo-  The  silver  plasma  frequency  is  13.8  x  1015  Hz  and  the  silver  plasma 

dephasing  is  0. 1075  x  1015  Hz.  The  energy  gap  Eg  of  the  quantum-dot  host  material  is  1 .927  eV 
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at  T  =  300  K.  Other  parameters,  including  T,  £sp,  Lx,  hjo,  Zo  and  fy ,  will  be  directly  indicated 
in  the  figures. 

Figure  3  presents  the  quantum  dot  absorption  coefficient  /3o(ft>sp)  for  an  SPP  field,  the  scat¬ 
tered  field  |Etot  —  Esp|  of  the  SPP  field,  and  the  energy-level  occupations  for  electrons  n,(  e  and 
holes  n.jth  with  £,  j  =  1,  2  as  functions  of  frequency  detuning  A7icosp  =  hcosp  —  (Eq  +  £i,e  +  £.\jt). 
A  dip  is  observed  at  resonance  A ti(Osp  =  0  in  the  upper-left  panel,  which  appears  to  become 
deeper  with  decreasing  amplitude  £sP  of  the  SPP  field  in  the  strong-coupling  regime  due  to  a 
decrease  in  the  saturated  absorption.  However,  this  dip  completely  disappears  when  Esp  drops 
to  25  kV/cm  in  the  weak-coupling  regime  due  to  the  suppression  of  the  photon-dressing  effect, 
which  is  accompanied  by  an  order  of  magnitude  increase  in  the  absorption-peak  strength.  The 
dip  in  the  upper-left  panel  corresponds  to  a  peak  in  the  scattered  field,  as  can  be  seen  from  the 
upper-right  panel  of  the  figure.  The  scattered  field  increases  with  frequency  detuning  away  from 
resonance,  corresponding  to  the  decreasing  absorption.  As  a  result,  a  minimum  appears  on  each 
side  of  the  resonance  in  the  scattered  field  in  the  strong-coupling  regime.  The  Maxwell-Bloch 
equations  couple  the  field  dynamics  outside  of  a  quantum  dot  with  the  electron  dynamics  inside 
the  dot.  At  £sp  =  125  kV/cm  in  the  lower-right  panel,  we  find  peaks  in  the  energy-level  occupa¬ 
tions  at  resonance,  which  are  broadened  by  the  finite  carrier  lifetime  as  well  as  the  optical  power 
of  the  SPP  field.  Moreover,  jumps  in  the  energy-level  occupations  can  be  seen  at  resonance  due 
to  Rabi  splitting  of  the  energy  levels  in  the  dressed  electron  states.  The  effect  of  resonant  phonon 
absorption  also  plays  a  significant  role  in  the  finite  value  of  «2,e  with  energy-level  separations 
£2,e  —  £i,e  ~  /A2().  However,  as  Tsp  decreases  to  25  kV/cm  in  the  lower-left  panel,  peaks  in  the 
energy-level  occupations  are  greatly  sharpened  and  negatively  shifted  due  to  the  suppression 
of  the  broadening  from  the  optical  power  and  the  excitonic  effect,  respectively.  Additionally, 
jumps  in  the  energy-level  occupations  become  invisible  because  the  Rabi-split  energy  gap  in 
this  case  is  much  smaller  than  the  energy-level  broadening  from  the  finite  lifetime  of  electrons 
(i.e.  severely  damped  Rabi  oscillations  between  the  first  electron  and  hole  levels). 

We  know  that  a  decrease  in  temperature  T  gives  rise  to  an  increase  in  the  crystal  bandgap  Eq. 
On  the  other  hand,  the  localization  of  an  SPP  field  (i.e.  an  exponential  decay  of  the  field  strength 
on  either  side  of  a  metallic  surface)  is  greatly  enhanced  when  the  SPP  frequency  ft)sp  approaches 
that  of  a  surface  plasmon.  As  a  result,  the  field  at  the  quantum  dot  is  expected  to  decrease 
as  T  is  reduced.  This  gives  rise  to  a  higher  absorption  coefficient  for  a  lower  temperature, 
as  shown  in  the  upper-left  panel  of  Fig.  4.  Interestingly,  although  the  suppressed  absorption 
coefficient  can  be  seen  from  /3o(ft)sp)  for  high  SPP-field  amplitudes,  as  shown  by  Eq.  (63), 
from  the  upper-right  panel  of  this  figure  we  find  the  resonant  peak  at  h(Osp  =  Eq  +  £i,e  +  £i  .h 
initially  increases  with  T  but  then  decreases  with  T  at  room  temperature.  This  subtle  difference 
demonstrates  the  effect  of  reduced  phonon  absorption  at  T  =  77  K  on  the  resonant  scattered 
field  by  the  factor  1  —  ne(t)  —  nh(f)  in  Eq.  (49).  Moreover,  the  strong  effect  of  the  suppressed 
optical-phonon  absorption  between  two  electron  energy  levels  at  77  K  is  clearly  demonstrated 
in  the  lower  panels  of  Fig.  4,  where  the  level  occupation  «2,e  becomes  negligible  at  T  =  77  K  in 
comparison  with  that  at  T  =  300  K. 

The  electron  thermal  dynamics  due  to  phonon  absorption  has  been  demonstrated  in  Fig.  4 
for  various  temperatures.  In  Fig.  5,  we  present  the  electron  dynamics  resulting  from  the  optical 
dephasing,  due  to  the  finite  lifetime  of  electrons,  at  different  energy-level  broadenings  Tiyo-  As 
Tr/o  is  increased  from  3  meV  to  7  meV,  the  dip  in  /3o(cosp)  at  resonance  is  suppressed,  leading  to 
a  single  peak  with  a  reduced  strength  and  an  increased  width,  as  shown  in  the  upper-left  panel 
of  the  figure.  This  increase  in  the  resonant  absorption  is  further  accompanied  by  an  enhanced 
resonant  peak  for  the  scattered  field  in  the  upper-right  panel  of  this  figure.  As  expected,  the 
energy-level  occupations  at  hjo  =  7  meV  become  much  broader  than  those  at  hyo  =  3  meV,  as 
displayed  in  the  lower  two  panels  of  the  figure. 
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Fig.  3.  Optical  absorption  coefficients  /3o  (cosp)  (upper-left  panel)  and  scattering  field  |Etot  — 
Esp|  at  the  quantum  dot  (upper-right  panel),  as  well  as  the  energy-level  occupations  for 
electrons  it(  e  and  holes  njh  (lower  panels)  with  l,  j  =  1,2,  as  functions  of  the  frequency 
detuning  AhdOsp  =  7ift)sp  —  (Eq  +£ijC  +  £ijt).  Here,  the  results  for  various  amplitudes  Esp  of 
an  SPP  field  with  frequency  are  presented  in  the  upper  panels,  along  with  a  comparison 
of  the  energy-level  occupations  for  Esp  =  25  and  125  kV/cm  in  the  lower  panels.  The  label 
xO.l  in  the  upper-left  panel  indicates  that  the  result  is  multiplied  by  a  factor  of  0.1. 


We  further  notice  that  the  effective  bandgap  Eq  +  £Pe  +  £|  .h  also  depends  on  the  size  Lx 
of  a  quantum  dot  due  to  the  quantization  effect,  and  the  effective  bandgap  will  increase  with 
decreasing  Lx.  The  size  effect  from  such  an  Lx  dependence  is  displayed  in  Fig.  6.  From  the 
upper-left  panel  of  Fig.  6,  we  find  that  the  peak  of  /3o(ft>Sp)  is  enhanced  as  Lx  is  reduced.  This 
phenomenon  is  connected  to  the  increased  localization  of  the  SPP  field  at  Lx  =  170  A  as  the 
SPP  frequency  approaches  the  saturation  part  of  its  dispersion.  Moreover,  the  dip  in  /3o(ft)Sp) 
is  lifted  somewhat  uniformly  at  the  same  time  due  to  a  decreased  n\(t)  from  the  enhanced 
Coulomb  and  phonon  scattering  at  Lx  =  170  A.  Here,  /3o(ft)sp)  is  proportional  to  the  population 
factor  1  —  n\(t)  —  n\(t),  as  can  be  seen  from  Eq.  (63).  Besides  the  slightly-reduced  resonant 
peak  strength  of  the  scattered  field  for  Lx  =  170  A  (also  resulting  from  the  enhanced  carrier 
scattering),  |Etot  —  Esp|  keeps  the  same  peak  position,  as  shown  in  the  upper-right  panel  of  the 
figure.  In  this  case,  |Etot  —  Esp|  at  the  dot  approaches  a  nonzero  value  at  resonance,  as  can  be 
seen  from  Eq.  (55),  and  tends  to  zero  rapidly  away  from  resonance.  Additionally,  «2,h  is  reduced 
forLt  =  170  A,  as  can  be  found  from  a  comparison  between  the  two  lower  panels  of  the  figure. 
This  is  attributed  to  the  reduced  phonon  absorption  between  two  hole  energy  levels. 

In  Figs.  4  and  6,  we  vary  the  localization  of  an  SPP  field  by  changing  the  effective  bandgap. 
Since  the  frequency  of  the  surface  plasmon  (saturated  dispersion  part)  is  proportional  to  the 
factor  of  \/\J\  +£d,  a  smaller  value  of  £(j  implies  a  higher  surface-plasmon  frequency  or  a 
reduced  localization  of  the  SPP  field.  We  verify  the  change  in  the  SPP  localization  by  observing 
the  upper  two  panels  of  Fig.  7,  where  the  absorption  peak,  as  well  as  the  resonant  scattered-field 
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Fig.  4.  /3o(fflsP)  (upper-left  panel)  and  |Etot  —  Esp|  (upper-right  panel),  as  well  as  ii(  e  and 
holes  rij  h  (lower  panels),  as  functions  of  ATitUsp.  Here,  the  results  for  three  different  temper¬ 
atures  T  =  300,  175  and  77  K  are  displayed  in  the  upper  panels,  along  with  a  comparison 
of  ri(  e  and  holes  njy t  for  T  =  300  and  77  K  in  the  lower  panels. 


peak,  become  much  stronger  as  £<i  is  increased  from  8  to  12  due  to  the  reduction  of  saturated 
absorption  for  a  lower  field  strength  at  the  quantum  dot.  Furthermore,  from  the  two  lower  panels 
of  this  figure  we  also  observe,  via  the  jumps  in  the  population  curves,  an  enhanced  Rabi-split 
energy  gap  in  the  electron  dressed  states  as  £d  is  reduced  from  12  to  10  due  to  the  enhanced 
field  strength  at  the  quantum  dot. 

In  the  presence  of  the  localization  of  an  SPP  field,  we  can  move  a  quantum  dot  closer  to 
a  metallic  surface  to  gain  a  higher  field  at  the  quantum  dot.  The  upper-left  panel  of  Fig.  8 
has  elucidated  this  fact,  in  which  a  larger  zo  corresponds  to  a  weaker  field,  and  then,  a  higher 
absorption  peak  due  to  the  reduction  of  saturated  absorption.  This  fact  is  also  reflected  in  the 
upper-right  panel  of  the  figure,  where  a  higher  resonant  scattered-field  peak  occurs  for  a  larger 
value  of  zo-  Atzo  =  510  A,  a  Rabi-split  energy  gap  at  resonance  is  clearly  visible  from  the  lower- 
left-panel  of  the  figure  for  electron  dressed  states.  Additionally,  at  zo  =  710  A,  by  entering  into 
a  weak-coupling  regime  for  a  weaker  field  at  the  dot,  we  find  sharpened  resonant  peaks  in  the 
energy-level  occupations  of  electrons  and  holes,  similar  to  the  observation  from  the  lower-left 
panel  of  Fig.  3. 

3.2.  Results  for  the  dressed  states  of  electrons 

In  the  second  part  of  the  numerical  calculations,  besides  the  parameters  given  in  the  first  sub¬ 
section,  we  have  fixed  Lx  =  210  A,  hyo  =  3meV,  zo  =  610  A  and  £<]  =  12.  Other  parameters, 
including  T,  Esp  and  A/icosp,  will  be  directly  indicated  in  the  figures.  Additionally,  A/icosp  is 
given  with  respect  to  the  energy  gap  at  T  =  300  K. 

From  the  left  panel  of  Fig.  9  we  find  a  strong  absorption  (positive)  peak  and  a  weak  gain  (neg- 
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Fig.  5.  /3o(®sp)  (upper-left  panel)  and  |Etot  —  Esp|  (upper-right  panel),  as  well  as  ri(  e  and 
holes  rij  h  (lower  panels),  as  functions  of  Ah(0Sp.  Flere,  the  results  with  different  energy- 
level  broadening  /i/o  =  3,  5  and  7  meV  are  shown  in  the  upper  panels,  along  with  a  com¬ 
parison  of  n(  e  and  holes  nj  ^  for  Tz/q  =  3  and  7  meV  in  the  lower  panels. 


ative)  peak  for  the  probe-field  absorption  coefficient  /3abs  ( 0)p )  due  to  a  quantum  coherence  effect 
from  the  electron  states  being  dressed  by  an  SPP  field.  In  the  strong-coupling  regime,  the  dis¬ 
persion  of  the  quantum-dot  e-h  plasmas  (dot-like  branch)  and  SPPs  (photon-like  branch)  form 
an  anticrossing  gap,  where  a  higher-energy  dot-like  branch  at  a  negative  frequency  detuning 
switches  to  a  photon-like  branch  for  a  positive  detuning.  The  positive  peak  is  associated  with  the 
absorption  of  a  probe-field  photon  by  a  quantum-dot  e-h  plasma,  while  the  negative  peak  relates 
to  the  process  with  absorption  of  two  photons  from  an  SPP  field  and  emission  of  one  probe- 
field  photon.  The  absorption  peak  is  significantly  reduced  by  saturation  at  Es p  =  lOOOkV/cm, 
and  the  gain  peak  is  suppressed  by  a  smaller  Rabi-coupling  frequency  at  Esp  =  250  kV/cm  (see 
the  inset  of  the  left  panel).  In  addition,  we  observe  from  the  right  panel  of  Fig.  9  that  two  Rabi- 
splitting-induced  side  emission  peaks  for  the  spontaneous  emission  Pem(co)  become  weaker 
and  closer  to  the  strong  central  peak  as  Esv  is  reduced  (see  the  inset  of  the  right  panel).  More¬ 
over,  the  strength  of  the  central  peak  due  to  the  coherent  conversion  of  an  absorbed  SPP-field 
photon  to  a  spontaneously-emitted  photon  (non-linear  optical  behavior)  is  slightly  reduced  at 
Esp  =  lOOOkV/cm  as  a  result  of  saturated  absorption  of  the  SPP  field. 

Figure  10  demonstrates  the  effect  of  frequency  detuning  Ahcosv  of  an  SPP  field  with  respect  to 
the  bandgap  of  a  quantum  dot.  The  switching  of  the  detuning  from  10  meV  to  — 10  meV  reveals 
the  corresponding  spectral-position  interchange  between  the  absorption  (dot-like  branch)  and 
the  gain  (photon-like  branch)  peaks  for  /3abs(&>p)  in  the  left  panel  of  the  figure.  The  Rabi  oscil¬ 
lations  between  the  first  electron  and  hole  energy  levels  are  weakened  with  increasing  |A/UOsp|. 
At  resonance  with  a  zero  detuning,  both  the  absorption  and  gain  peaks  are  suppressed  by  very 
strong  Rabi  oscillations.  This  detuning  also  shifts  the  emission  peaks  correspondingly  because 
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Fig.  6.  /3o  ( CDsp )  (upper-left  panel)  and  |Etot  —  Esp|  (upper-right  panel),  as  well  as  ng  e  and 
holes  tijj,  (lower  panels),  as  functions  of  AhcoSp.  Here,  the  results  for  three  different  sizes 
Lx  =  210,  190  and  170A  of  a  quantum  dot  are  shown  in  the  upper  panels,  along  with  a 
comparison  of  nge  and  holes  nj for  Lx  =  210  and  170  A  in  the  lower  panels. 


of  the  coherent  conversion  of  an  SPP-field  photon  to  a  spontaneously-emitted  one,  as  can  been 
seen  from  the  right  panel  of  this  figure.  Moreover,  the  central  peak  is  weakened  and  the  two 
side  peaks  are  enlarged  at  resonance  as  a  result  of  energy  transfer  to  the  side  peaks  by  strong 
coupling  and  enhanced  Rabi  oscillations,  respectively. 

Since  the  temperature  affects  the  crystal  bandgap  energy  Eq,  by  changing  the  temperature 
we  are  able  to  scan  the  detuning  A /ift)sp  of  the  SPP  field  with  a  fixed  SPP  frequency  /zO)sp 
from  negative  to  positive  or  vice  versa.  This  leads  to  a  spectral-position  interchange  between 
the  absorption  and  gain  peaks,  similar  to  Fig.  10.  The  results  in  Fig.  1 1  prove  such  an  expected 
feature  by  increasing  T  from  250  to  300  K  in  steps  of  5  K.  Technically,  changing  the  temperature 
in  the  experiment  is  much  easier  than  changing  the  tuning  of  a  laser  frequency  over  a  large 
range.  Here,  the  shift  of  the  central  peak  in  the  right  panel  of  the  figure  directly  reflects  the 
variation  of  the  SPP-field  detuning  with  T.  Furthermore,  the  interchange  between  the  dot-like 
and  photon-like  modes  with  T  in  the  left  panel  can  be  regarded  as  direct  evidence  for  the 
existence  of  an  anticrossing  energy  gap  resulting  from  a  strongly-coupled  e-h  plasma  and  SPP 
field  or  coupled  e-h  plasmas  and  surface  plasmons. 

3.3.  Time-resolved  optical  spectra 

In  our  previously  presented  numerical  results,  we  only  showed  steady-state  dynamics  of  photo- 
excited  e-h  plasmas  in  a  quantum  dot  by  using  a  continuous  SPP  field,  where  the  effects  of 
both  phonon  scattering  and  e-h  pair  radiative  recombination  are  combined  with  each  other. 
Using  a  laser  pulse  to  launch  a  pulsed  SPP  field,  we  are  able  to  study  the  dynamics  of  phonon 
scattering  (narrow  pulse)  as  well  as  the  dynamics  of  e-h  pair  radiative  recombination  (wide 
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Fig.  7.  /3o(fflsP)  (upper-left  panel)  and  |Etot  —  Esp|  (upper-right  panel),  as  well  as  ng  e  and 
holes  tij  h  (lower  panels),  as  functions  of  Afuusp.  Flere,  the  results  for  three  dielectric  con¬ 
stants  of  a  cladding  layer,  with  Ej  =  8,  10  and  12,  are  displayed  in  the  upper  panels,  along 
with  a  comparison  of  n(  e  and  holes  nj for  =  12  and  10  in  the  lower  panels. 


pulse),  separately.  Dynamically,  phonon  scattering  becomes  effective  only  after  a  characteristic 
time  (around  1  ps),  its  effect  can  be  seen  from  a  significant  increase  of  «2,e  in  our  system.  Figure 
12  displays  the  results  for  /3o(coSp)  (upper-left),  |Etot  —  Esp|  (upper-right),  n\,e  (lower-left)  and 
«2.e  (lower-right)  for  various  detection  times  To  in  the  presence  of  a  narrow  laser  pulse  (with 
pulse  width  Tp  =  500  fs  and  peak  value  Esp  =  500kV/cm)  which  is  turned  on  at  t  =  0.  We 
see  from  Fig.  12  that  /3o(ft>Sp)  starts  with  a  dip  for  the  dressed  state  at  resonance,  then  shifts  to  a 
single  peak  (at  half-pulse  width)  due  to  a  suppression  of  the  photon-dressing  effect.  It  eventually 
becomes  a  single  peak  plus  a  shifted  dip  after  the  pulse  has  passed  due  to  formation  of  resonant 
peaks  in  «i.e  and  n\^.  Correspondingly,  |Etot  —  Esp|  starts  by  showing  a  non-resonant  behavior 
with  a  relatively  large  magnitude,  then  shifts  to  a  quasi-resonant  behavior,  and  finally  looks  like 
suppressed  resonant  behavior  with  a  peak  at  and  dips  on  both  sides  of  A h(Osp  =  0.  The  resonant 
build  up  of  nj  e  after  To  >  500  fs  can  also  be  verified  from  this  figure,  which  is  accompanied  by 
the  start  of  significant  phonon  absorption  after  To  >  1  ps. 

Technically,  detecting  dynamics  of  photo-excited  e-h  plasmas  by  using  another  time-delayed 
weak  probe  field  is  much  more  feasible,  as  shown  in  Fig.  13.  From  the  left  panel  of  this  figure, 
we  find  that  /3abs(®p)  starts  with  a  pair  of  positive  absorption  and  negative  gain  peaks  due  to 
a  very  strong  photon  dressing  effect  for  the  delay  times  l,i  =  60  and  120  fs.  This  is  changed 
to  a  strong  absorption  peak  plus  a  very  weak  gain  peak  at  T cj  =  240  fs.  At  the  end,  /3ahs(cop) 
becomes  independent  of  Tj,  indicating  that  a  linear  optical-response  regime  has  been  reached. 
On  the  other  hand,  from  the  right  panel  of  this  figure,  we  see  that  the  central  peak  of  Pem(o)  ')  is 
gradually  built  up  with  increasing  Tj  due  to  enhanced  n\e  and  n  |  h  around  resonance,  while  two 
side  peaks  become  weakened  and  disappear  at  the  same  time  due  to  weakened  Rabi  oscillations. 
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Fig.  8.  /3()(ffl sP)  (upper-left  panel)  and  |Etot  — Esp|  (upper-right  panel),  as  well  as  nge  and 
holes  n  jj,  (lower  panels),  as  functions  of  A/zft»sP.  Flere,  the  results  for  different  vertical 
distances  of  a  quantum  dot  from  a  metallic  surface  are  displayed  in  the  upper  panels  for 
zo  =  510,  610  and  710  A,  along  with  a  comparison  of  nge  and  holes  nj  j,  for  zo  =  510  and 
710  A  in  the  lower  panels. 


H'(EG+E1,e+£1,h)  (meV) 


h(0  -  (EG+e1  e+s1  h)  (meV) 


Fig.  9.  Probe-field  absorption  coefficient  /3abs(ft>p)  (left  panel)  and  spontaneous  emission 
of  a  quantum  dot  Pem(®)  (right  panel)  as  functions  of  7i(0p  —  (Eg  +  £| ,e  +  £i,h)  and  h(0  — 
(Eq  -f  £i  e  +  £i,h),  respectively,  are  presented.  Here,  comparisons  of  the  results  with  three 
values  of  SPP-field  amplitudes  Esp  =  250,  500  and  lOOOkV/cm  are  given.  The  insets  of 
both  panels  are  enlarged  views  of  the  peaks. 
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Fig.  10.  /3abs(ft*p)  (left  panel)  and  Pem(ft>)  (right  panel)  as  functions  of  TicOp  —  (Eg  +  £i,e  + 
£ph)  and  Tico—  (Eg  +  £i,e  +  £i,h)>  respectively,  are  displayed.  Flere,  we  show  the  compar¬ 
isons  with  three  values  of  SPP-field  detunings  Ah(Osv  =  0  and  ±10meV.  The  label  x  100 
for  zero  SPP-field  detuning  indicates  the  result  is  multiplied  by  a  factor  of  100. 


Fig.  11.  /3abs(Wp)  (left  panel)  and  Pem(ft>)  (right  panel)  as  functions  of  h(Op  —  (Eg  +  £i.e  + 
£yh)  and  tico  —  (Eg  +  £i.e  +  £i.h)>  respectively,  are  shown.  Flere,  a  series  of  offset  curves 
are  displayed  for  various  temperatures  from  T  =  250  K  (bottom  black  curves)  to  T  =  300  K 
(top  blue  curves)  in  steps  of  5  K. 


Interestingly,  we  also  find  that  the  central  peak  of  Pem(co)  slightly  decreases  at  TtJ  =  1  ps,  which 
agrees  with  the  observed  start  of  significant  phonon  absorption  seen  in  the  lower-left  panel  of 
Fig.  12. 

In  order  to  explore  the  dynamics  of  e-h  pair  radiative  recombination  in  our  system,  a  wide 
pulse  with  a  full -pulse  width  around  300 ps  is  required,  as  displayed  in  Fig.  14.  From  the  upper- 
middle  panel  of  this  figure,  we  find  that  /3o(fflsp)  starts  with  a  resonant  dip  due  to  a  strong  photon 
dressing  effect,  then  shifts  to  a  sole  peak  at  A/i ft)sp  =  0  as  To  >  400  ps  where  a  steady  state  is 
almost  reached  in  the  linear-response  regime.  Accordingly,  the  level  populations  n\  e  and  «2,e  in 
the  lower  two  panels  show  a  transition  from  an  initial  non-resonant  behavior  to  a  final  resonant 
behavior.  This  is  accompanied  by  dramatically  reduced  level  populations  due  to  the  start  of  a 
radiative  recombination  process  for  photo-excited  e-h  pairs. 

Recombination  dynamics  for  e-h  plasmas  can  also  be  demonstrated  clearly  by  the  time- 
delayed  probe-field  absorption  as  well  as  by  the  time-resolved  spontaneous  emission,  as  shown 
in  Fig.  15.  As  presented  in  the  left  panel  of  this  figure,  we  find  that  the  initial  weak  absorption 
and  gain  peaks  (see  the  inset)  in  j3abs(ft>p)  occur  at  xd  =  200  ps  and  are  replaced  by  a  strong 
single  absorption  peak  due  to  a  suppressed  photon  dressing  effect  and  phase-space  blocking. 
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Fig.  12.  j3o(ft)sp)  (upper-left  panel)  and  |Etot  —  Esp|  (upper-right  panel),  as  well  as  n\  e 
and  holes  «2  e  (lower  panels),  as  functions  of  A7iC0sP  at  different  detection  times  To  for  a 
Gaussian-shape  laser  pulse  with  pulse  width  Tp  =  500  fs.  Flere,  £sp  =  500kV/cm  is  taken, 
and  the  other  parameters  are  the  same  as  those  in  Fig.  3.  The  labels  x  10  and  x  105  in  the 
upper-right  panel  indicate  that  the  results  are  multiplied  by  factors  of  10  and  10s,  respec¬ 
tively. 


Fig.  13.  /3abs(ft*p)  (left  panel)  and  Pe m(ft>)  (right  panel)  as  functions  of  TicOp  —  (Eq  +  £i  e  + 
£j  h)  and  ti(0  —  (Eq  +  £i  e  +  E\  h),  respectively,  at  different  delay  times  T^  for  a  Gaussian- 
shape  laser  pulse  with  pulse  width  Tp  =  500  fs.  Flere,  Esv  =  500  kV/cm  and  A/iftijp  =  5  meV 
are  chosen,  and  the  other  parameters  are  the  same  as  those  in  Fig.  3. 
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Fig.  14.  /3o  ( coSp )  (upper-middle  panel),  n\  e  and  holes  «2  e  (lower  two  panels),  as  functions 
of  Ah(0Sp  at  different  detection  times  To  for  a  wide  laser  pulse  with  pulse  width  Tp  =  300  ps. 
Here,  Esp  =  500kV/cm  is  taken,  and  the  other  parameters  are  the  same  as  those  in  Fig.  3. 
The  labels  x0.05  in  the  upper-middle  panel  indicate  that  the  results  in  the  upper  panel  for 
To  =  60  and  120 ps  are  multiplied  by  a  factor  of  0.05. 


On  the  other  hand,  from  the  right  panel  of  the  same  figure,  we  see  that  the  initial  central  peak 
in  Pe m(ft))  is  increased  very  rapidly  due  to  accumulation  of  photo-excited  e-h  pairs  and  ac¬ 
companied  by  the  reduction  of  two  side  peaks  resulting  from  the  weakened  Rabi  oscillations. 
Importantly,  the  very-strong  central  peak  in  Pem(co)  is  significantly  reduced  at  T,/  =  200 ps, 
indicating  the  start  of  a  radiative-recombination  process  for  photo-excited  e-h  plasmas.  This  re¬ 
combination  process  is  continuously  enhanced  with  the  increasing  delay  time  xd  and  suppresses 
the  central  peak  in  Pem((o)  after  t d  >  400 ps  due  to  draining  out  the  photo-generated  electrons 
and  holes  at  the  same  time. 

4.  Conclusions  and  remarks 

In  conclusion,  we  have  demonstrated  the  possibility  of  using  a  SPP  field  to  control  the  optical 
gain  and  absorption  of  another  passing  light  beam  due  to  their  strong  nonlinear  field  coupling 
mediated  by  electrons  in  the  quantum  dot.  We  have  also  predicted  the  coherent  conversion  of  a 
surface-plasmon-field  photon  to  a  spontaneously-emitted  free-space  photon,  which  is  simulta¬ 
neously  accompanied  by  another  pair  of  blue-  and  red-shifted  photons. 

Although  we  studied  only  the  coupling  of  a  SPP  field  to  a  single  quantum  dot  in  this  paper  for 
the  simplest  case,  our  formalism  can  be  generalized  easily  to  include  many  quantum  dots.  The 
numerically-demonstrated  unique  control  of  the  effective  photon-photon  coupling  by  the  quan¬ 
tum  dot  can  be  used  for  constructing  an  optical  transistor,  where  the  ‘gate’  photon  controls  the 
intensity  of  its  ‘source’  light  beam.  These  optical  transistors  are  very  useful  for  speeding  up  and 
improving  the  performance  of  fiber-optic  communication  networks,  as  well  as  for  constructing 
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Fig.  15.  /3abs(fttp)  (left  panel)  and  Pe m(co)  (right  panel)  as  functions  of  ftcdp  —  (Eg  +  £i.e  + 
£i,h)  and  hco  —  (Eg  +  £iie  +  £i  h).  respectively,  at  different  delay  times  T</  for  a  Gaussian- 
shape  laser  pulse  with  pulse  width  Tp  =  300  ps.  Here,  £sp  =  500  kV/cm  and  A/7ffi>sp  =  5  meV 
are  assumed,  and  the  other  parameters  are  the  same  as  those  in  Fig.  3.  The  inset  in  the  upper- 
left  panel  shows  an  enlarged  view  for  the  gain  and  absorption  peaks  for  small  delay  times. 
The  labels  x  10  for  =  50  and  lOOps  in  the  left  panel,  as  well  as  the  label  x0.02  in  the 
right  panel,  indicate  that  the  results  are  multiplied  by  a  factor  of  10  and  a  factor  of  0.02, 
respectively. 


quantum  information  and  developing  optical  digital  computers. 

Furthermore,  instead  of  a  resonant  coupling  to  the  lowest  pair  of  electron-hole  energy  lev¬ 
els,  we  may  select  the  surface-plasmon  frequency  for  resonant  coupling  to  the  higher  pair  of 
electron-hole  levels.  In  such  a  case,  the  optical  pumping  from  the  intense  surface-plasmon  near- 
field  could  create  a  population  inversion  with  respect  to  the  ground  pair  of  electron-hole  levels 
by  emitting  thermal  phonons,  leading  to  a  possible  lasing  action  if  the  optical  gain  can  over¬ 
come  the  metal  loss  for  the  surface  plasmons.  Such  a  surface-plasmon  quantum-dot  laser  would 
have  a  beam  size  as  small  as  a  few  nanometers  beyond  the  optical  diffraction  limit,  and  it  is  ex¬ 
pected  to  be  very  useful  for  spatially-selective  illumination  of  individual  molecules  or  neuron 
cells  in  low-temperature  photo-excited  chemical  reactions  or  optogenetics  and  neuroscience. 

A.  Results  Related  to  the  Electronic  States  of  a  Quantum  Dot 

We  have  employed  a  box-type  potential  with  hard  walls  to  model  a  quantum  dot,  which  is  given 
by 


j°,  0p,<L,for/=l,2,3 

[  °° ,  otherwise 

where  the  position  vector  r  =  (x\,X2,x-i),  and  L\ ,  Ln  and  Lj,  are  the  widths  of  the  potential  in 
the  xi,X2  and  *3  directions,  respectively.  The  Schrodinger  equation  for  a  single  electron  or  hole 
in  a  quantum  dot  is  written  as 


2m* 


d2  d 2  d2  . 

[dx2  +  dx2  +  dx2+V^) 


V(r)  = 


where  the  effective  mass  m*  is  mL 
associated  with  Eq.  (74)  is  found  to  be 


for  electrons  or  rn'b  for  holes. 


e  yr(r)  ,  (74) 

The  eigenstate  wave-function 
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Wnltn2,n3  (r) 


(75) 


which  is  same  for  both  electrons  and  holes,  and  the  eigenstate  energy  associated  with  Eq.  (74) 
is 


£n\ 


2m* 


(76) 


where  the  quantum  numbers  «i ,  n.2,  «3  =  1,  2,  •  •  •. 

By  using  the  calculated  bare  energy  levels  in  Eq.  (76),  the  dressed  electron  (A^)  and  hole 
(A£)  energy  levels  under  the  rotating  wave  approximation  take  the  forms  [2] 


Xea{(o\t)=X^{co\t)=  < 


3  ( h a  +  \J[£g{T)  +  £„  +  £„  —  h(o}2  +  4|-#®hcf(t)|2^ 

if  h(0  <  (T1 )  H-  H- 

3  —  \J{^g{T)  +  ea  +  ea  ~  ^®]2  +4|^Jha(t)|2 

if  /ift)  >  (^(r)  +£«  +£« 


(77) 


where  the  composite  index  a  =  {n i,  ni.  n{\.  Moreover,  we  obtain  the  energy  levels  of  dressed 
electrons  esa((o\t)  =  A^( to |f)  +  (e„  —  eJj)/2  and  e|(co|r)  =  e|  +  <%(r)/2  for  £  ^  a.  Similarly, 
we  obtain  the  energy  levels  of  dressed  holes  e^(to|f)  =  A^(to|f)  +  (ejj  —  e®)/2  and  Ej(co\t)  = 
E/1  +  <£g(T)/2  for  j  ±  a. 

Based  on  the  calculated  wave-functions  in  Eq.  (75),  the  form  factors  introduced  in  Eqs.  (11) 
and  (12)  can  be  obtained  from 


^ni,n2,ny,n!vn'2,n'3^  ~  ^ni  ,n2lny,n[,n'2,n'3^  >  ^8) 


(1) 


a(2) 


9(3) 


where  the  wave  vector  q  =  (<71 ,  <72 ,  <?3 )  and  we  have  introduced  the  following  notation  for  j  = 
1,2,3 


Kj^j)=[j^)  JdxjfPm  sin 


0 


sin 


n)K 


(79) 


Moreover,the  overlap  of  the  electron  and  hole  wave-functions  in  this  model  can  be  easily  cal¬ 
culated  as 


/ 


^r<,„2,„3( r)  ¥^n'2y3 (r)  =  5ni,„'  K2J2  S„3,„' 


(80) 


The  interband  dipole  moment  dc  v  =  dcve(i  at  the  isotropic  T-point,  which  is  defined  in  Eq.  (32), 
can  be  calculated  according  to  the  Kane  approximation  [35, 36] 


l/c,V  — 


e2Tr 


2mQ<gG(T)  \ni 


m  _  1 


(81) 


Furthermore,  the  direction  of  the  dipole  moment  e<j  is  determined  by  the  quantum-dot  energy 
levels  in  resonance  with  the  photon  energy  Tift). 
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